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Abstract. Highly oriented solid-supported lipid membranes in stacks of controlled number N ~ 16 (oligo- 
membranes) have been prepared by spin-coating using the uncharged lipid model system 1,2-dimyristoyl- 
sn-glycero-3-phosphochohne (DMPC). The samples have been immersed in aqueous polymer solutions 
for control of osmotic pressure and have been studied by X-ray reflectivity. The bilayer structure and 
fluctuations have been determined by modelling the data over the full g-range. Thermal fluctuations are 
described using the continuous smectic hamiltonian with the appropriate boundary conditions at the 
substrate and at the free surface of the stack. The resulting fluctuation amplitudes and the pressure- 
distance relation are discussed in view of the inter-bilayer potential. 

PACS. 61.10.Kw X-ray reflectometry (surfaces, interfaces, Aims) - 87.16.Dg Membranes, bilayers, and 
vesicles - 87.15.Ya Fluctuations 


1 Introduction 

A quantitative understanding of the structure, fluctua¬ 
tions, interaction potential and elasticity properties of lipid 
membranes, which represent model systems for biological 
membranes, has been the goal of many theoretical and 
experimental studies. Theoretically, they have been stud¬ 
ied as paradigmatic examples of quasi two-dimensional 
macromolecular structures governed by bending rigidity 
[T]. In aqueous solution, lipid bilayers assemble into stacks 
governed by distinct inter-bilayer interactions. A number 
of seminal studies using high resolution synchrotron X- 
rays have been published on these systems [11[21I11[S]. In 
these studies, isotropic aqueous dispersions of multilamel- 
lar vesicles have been studied as a function of temperature 
T and/or osmotic pressure U. Detailed quantitative in¬ 
formation on the interaction potentials and the elasticity 
properties has thus been derived, see e.g. mm- How¬ 
ever, information is lost in these small-angle scattering 
studies due to crystallographic powder averaging. In the 
analysis, assumptions must therefore be made on the na¬ 
ture of the correlation functions in the framework of linear 
smectic elasticity theory, leading to the Caille model [6] 
and related theories, see for instance [7]. In order to over¬ 
come the limitations of powder averaging, it is desirable 
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to work with aligned systems of lipid bilayers 
[H]. Under the same conditions of temperature and hy¬ 
dration, thermal fluctuations are not as strong for aligned 
systems as in bulk studies due to the boundary condition 
at the flat substrate, enabling a higher resolution in p{z). 
It is also advantageous to fit the data continuously over a 
large range of momentum transfer q, as has been shown 
for isotropic solutions and for oriented stacks and 
not only in the vicinity of the Bragg peaks arising from 
the multilamellar structure. However, for aligned systems 
of multilamellar membranes, satisfactory fits of the entire 
reflectivity curves and the formulation of a proper statis¬ 
tical model as well as of a scattering theory are still quite 
difficult. Note that the best published fits are for systems 
consisting of monolayers at the air-water interface, for sin¬ 
gle bilayers or for a free-floating bilayer system developed 
by Fragneto and coworkers m- In multilamellar systems 
on the other hand, the reflectivity signal is typically much 
more complex and structured. As we show here, structural 
parameters of the bilayer, interaction and fluctuation pa¬ 
rameters can be deduced from these curves, and are com¬ 
pared with the literature. 


Most previous studies on aligned multilamellar mem¬ 
branes suffered from a lack of control concerning the num¬ 
ber of bilayers N, the sample homogeneity and also the bi¬ 
layer hydration. Building on recent progress in the prepa¬ 
ration of bilayers on solid support by spin-coating El, 
we use so-called oligo-membranes with a reduced num¬ 
ber of bilayers N ~ 10 — 20 resulting in very structured 
and well resolved reflectivity signals. We have developed 
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a model for the thermal fluctuations calculated for the 
proper boundary conditions (rigid substrate and free up¬ 
per surface), which gives the fluctuation amplitudes 
of the bilayers as a function of their position in the stack, 
n = 1,7V. The values for (m^) are then inserted in the mul- 
tilamellar structure factor, along with a decreasing cover¬ 
age function (see below). The density profile p{z) is pa¬ 
rameterized in terms of its Fourier coefficients |15j . This 
approach gives for the first time an agreement with the 
measured reflectivity curve of multilamellar membranes 
over the full range of Qz up to typically qz — 0.7 A and 
over about seven orders of magnitude in the measured sig¬ 
nal. 

The main experimental parameter in this study is the 
osmotic pressure. The classical osmotic stress (OS) tech¬ 
nique as developed by Parsegian and coworkers [16] is 
widely used for the measurement of force-distance curves 
in colloidal systems. The osmotic pressure imposed to a 
lamellar phase controls the interaction force experienced 
by the membranes across the water layer by setting the 
chemical potential of the water molecules in the inter¬ 
bilayer solution. Pressure-distance relations can be easily 
determined, e.g. if the lamellar periodicity d is measured 
by X-ray scattering at different pressures U. 

In this study we use a variant of the OS technique 
where the oriented bilayers are put in direct contact with 
the osmotic stress solution [18]. We have verified that 
the high molecular weight polymer with a radius of gyra¬ 
tion larger than d does not penetrate the lamellar phase. 
We emphasize that the osmotic pressure is one of the 
most important parameters in biomolecular systems, since 
biomolecular assemblies in the cell are mostly exposed to 
varying U while T is often constant. Therefore, it is of 
great importance to study bilayer structure and elastic 
properties such as bilayer bending rigidity k or bilayer- 
bilayer interaction parameters as a function of osmotic 
pressure 11. 

From the relation between the osmotic pressure U and 
the lamellar spacing d the inter-bilayer interaction poten¬ 
tials can be determined. It is generally accepted that in 
charge neutral systems two main molecular interaction 
forces are dominant, in addition to the effective attrac¬ 
tive interaction by osmotic pressure : a repulsive hydra¬ 
tion potential /hyd(dm) and the attractive van der Waals 
potential /vdw(du,) so that the total interaction in J/m^ 
is given by f{du,) = fhyd + /vdw +11 d^i, defining the equi¬ 
librium distance (water layer thickness) d^ as illustrated 
in Fig. |I]B). As discussed below, it is important to take 
the correct form for /vdw as derived in HZ], without the 
conventional half-space approximation. It can be argued 
that steric (Helfrich) repulsion forces have to be added 
to the molecular forces in a mean field approach. Here, 
however, we will assume that thermal fluctuations in thin 
films of relatively stiff phospholipids do not have a signif¬ 
icant effect on the inter-bilayer interactions, in particular 
since the flat boundary suppresses long range fluctuations 
in the film. The paper is organized as follows : After this 
introduction, section [5] presents some experimental details 
while the statistical model and data analysis are presented 



(B) 



Fig. 1. A): Sketch of the experimental setup where a solid- 
supported membrane stack is hydrated in an aqueous polymer 
solution. Arrows indicate the incoming and scattered X-ray 
beams. B): Interaction potential with van der Waals contri¬ 
bution according to Fenzl [IT] (solid line) and with general 
approximation (dashed). The (dotted) parabola illustrates the 
harmonic approximation to the potential which enters in the 
smectic elasticity theory. The arrow shows a typical fluctuation 
amplitude of a membrane inside a 16 bilayer stack. 


in section 13] Section [J] presents the results, followed by a 
section on the interaction potentials and the conclusions 
in section i] 


2 Materials and Methods 

2.1 Samples and Environment 

Highly oriented oligo-membranes were prepared using the 
spin-coating method M- The uncharged lipid 1,2-dimy- 
ristoyl-sn-glycero-3-phosphocholine (DMPC) was bought 
from Avanti (Alabaster, AL, USA) and used without fur¬ 
ther purification. The lipid was dissolved in chloroform 
at a concentration of 10 mg/ml. An amount of 100 ^1 
of the solution was pipetted onto carefully cleaned sili¬ 
con substrates of a size of 15 x 25 mm^ cut from stan¬ 
dard commercial silicon wafers. The substrate was then 
immediately accelerated to rotation (3000 rpm), using a 
spin-coater. After 30 seconds the samples were dry and 
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subsequently exposed to high vacuum to remove any re¬ 
maining traces of solvent. The samples were then stored at 
4°C until the measurement. For the X-ray measurements 
the samples were hydrated in a stainless steel chamber 
m with kapton windows, which can be filled with wa¬ 
ter or with polymer solution to control the level of hy¬ 
dration by osmotic pressure, see the sketch in Fig. [TJ\. 
Temperature was controlled by an additional outer cham¬ 
ber at r = 40°C. The polymer polyethyleneglycol (PEG) 
of molar weight 20000 Da was bought from Fluka and 
used without further purification. PEG was dissolved in 
ultrapure water (Millipore, Billerica, Mass.) at the con¬ 
centrations 1.5 %, 2.9 %, 3.6 %, 5.8 %, 9 %, 12.1 %, 
14.2 % and 25 % (wt. percent). The corresponding os¬ 
motic pressure values were taken from the literature. The 
data was obtained from the web site of the Membrane 
Biophysics Laboratory at the Brock University in Ganada 
(http://aqueous.labs.brocku.ca/osfile.html). The value for 
the osmotic pressure of PEG 20000 solutions is only avail¬ 
able at 20 °C. At 40 °G, temperature at which the experi¬ 
ments were performed, the pressure can be expected to be 
somewhat lower. However, the temperature coefficient is 
small [MBS], and the corresponding discrepancy smaller 
than the error bars in Fig below, which is used in the 
analysis of the interaction forces. 

An important issue in our method of direct contact of 
the lamellar phase with the polymer solution is that of 
possible interpenetration of the multilamellar phase. Ex¬ 
periments on polymer containing lyotropic lamellar phases 
have shown that polymers can enter the water layer 
in between charged bilayers even if the radius of gyration 
is larger than the water layer. However, in the present case 
of neutral polymer in neutral bilayers, the amount of poly¬ 
mer in between the lamellae is negligible. The experimen¬ 
tal proof is given by (a) the density profile, which shows 
no deviation from the water density in between the bilay¬ 
ers, and (b) the force distance curve itself which shows no 
indication of such an effect. 


2.2 X-ray experiment 



Fig. 2. Specular reflectivity scan (open symbols) and offset 
scan (red solid curve) for a sample immersed in pure water. 
The "true specular” contribution (see Appendix) is given by 
the difference of the two curves. In dotted line we also show the 
Fresnel reflectivity profile corresponding to the same critical 
angle qc- Inset : Rocking scan at the position of the second 
Bragg peak. The arrow shows the position of the offset scan. 


scan on the second Bragg peak, illustrating the separation 
between ’’true specular” and diffuse components. 


3 Model and data analysis 

3.1 Reflectivity 


In the semi-kinematic approximation the reflectivity of 
a structured interface can be expressed by the so-called 
master-formula of reflectivity [H] as : 


R{qz) 


Rpiqz) ■ 


1 

Pl2 



dz 


2 


( 1 ) 


The X-ray reflectivity measurements presented here were 
carried out at the bending magnet beamline D4 of HA- 
SYLAB/DESY in Hamburg, Germany. At D4, a single¬ 
reflection Si(lll) monochromator was used to select a 
photon energy of 19.92 keV, after passing a Rh mirror to 
reduce higher harmonics. The chamber was mounted on 
the z-axis diffractometer, and the reflected beam was mea¬ 
sured by a fast scintillation counter (Cyberstar, Oxford), 
using computer-controlled aluminum absorbers which at¬ 
tenuate the beam at small Qz to prevent detector satura¬ 
tion. Incident and exit beams were defined by a system of 
several motorized slits. The data is corrected for decreas¬ 
ing electron ring current and the diffuse contribution by 
subtraction of an offset scan. Finally, an illumination cor¬ 
rection is performed. A typical measurement (reflectivity 
and offset scans) is shown in Figured along with the cor¬ 
responding Fresnel reflectivity. The inset shows a rocking 


where Rp denotes the Fresnel reflectivity of the sharp in¬ 
terface and p{z) is the intrinsic electron density profile, 
whereas pi 2 is the total step in electron density between 
the two adjoining media. The electron density profile of 
a solid-supported oligo-membrane stack, consisting of N 
membranes in water, can then be written as : 


Piz) = ^ erfc 


N-1 


+ > po{z-nd + Un), (2) 


n—0 


with erfc(z) being the complementary error function and 
a the rms substrate roughness. pi 2 = psi ~ PH 2 O is the 
difference in electron density between the substrate and 
water, do is the distance between the substrate and the 
midpoint of the first bilayer and Pq{z) is the electron den¬ 
sity profile of one bilayer in the stack. Thermal membrane 
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fluctuations are considered in terms of the displacement 
function = u{v\\,z = nd) of the position of the n-th 
membrane from its average position z = nd along z. Re¬ 
placing the electron density profile m into o and taking 
the ensemble average yields : 


R{qz) = Rpiqz) 

2 2 -^”1 / 

- 2i ■ e~^ ^ f{n) (• sin(g^(do + nd)) 


n—0 

N-1 


m,n—0 


The first summand represents the reflectivity of the sub¬ 
strate. The second is a cross-term and represents inter¬ 
ference effects between the substrate and the membrane 
stack. The third summand is the product of the form fac¬ 
tor \Ff{qz)\^ containing the structural information about 
one bilayer in the stack, and the structure factor, repre¬ 
senting the periodic structure in the stack. The fluctua¬ 
tions are described by the correlation function (unUm)- In 
specular reflectivity only the self-correlation function (u^) 
is important (see Appendix). 


3.2 Correlation Function 


The self-correlation function of the membrane fluctuations 
can be calculated from linear smectic elasticity theory 
based on a continuous model m- The complete theory 
and calculations are described in [53]. Here only the es¬ 
sentials, which are important for specular reflectivity shall 
be given. The linearized free energy can be written as a 
function of the displacement M(r||, z) as : 



/ du{rn,z) 
dz 


2 

-f A:(Z\||u(r||,z))^ 


(3) 


with K = n/d the bending modulus and B the compres¬ 
sion modulus in the stack. We neglect the surface tension 
between the lipid stack and the solvent. The discrete struc¬ 
ture of the stack consisting of N bilayers is taken into ac¬ 
count by expanding u over N independent modes. Also, we 
are only interested in the fluctuation amplitude at the po¬ 
sition of the bilayer midpoints, denoted by = u{z = nd) 
for the n-th bilayer. From the equipartition theorem one 
can calculate the correlation function of the membrane 
fluctuations (u^), which reads : 


(ul) =v{- 

TT 


2 N 

^2,-1 


■ sin 


2 j - 1 n 

-TT- 

2 N 


(4) 


with the conventional Caille factor rj = Figure 

Hv shows the function for a 16 membrane stack with a 
typical 77 value for DMPC membranes at partial hydra¬ 
tion. In contrast to oligo-membranes which are partially 


hydrated from water vapor, oligo-membranes immersed in 
excess water or polymer solution exhibit defects which re¬ 
sult in decreasing layer coverage with increasing distance 
from the substrate. This effect is evidenced experimen¬ 
tally by the suppression of thickness oscillations (Kiessig 
fringes) in the reflectivity curves. Thickness oscillations 
are typically observed in vapor-hydrated samples but are 
significantly reduced or suppressed in samples immersed in 
aqueous solution [m[55] . In the model we take this effect 
into account by multiplying the contribution of each mem¬ 
brane in the structure-factor with an empirical coverage- 
factor /(n) = [1 — (n/Af)“]^, where a parameterizes the 
decaying density due to decreasing coverage. A typical 
experimental value is a ~ 1.7. Note that both a and the 
number of bilayers N are ht parameters. 

Since the form factor \Ff{qz)\^ of the membranes con¬ 
sists of the squared Fourier transform of the z-derivative 
of the electron density profile of the membrane, it is conve¬ 
nient to express the profile in terms of normalized Fourier 
coefficients [15] 


No r) 

P{z) ='^ fm- COS ■ P 12 +P- Pw, 

m—l 

with p being the average electron density of the membrane 
stack and pw the electron density of water. For DMPC, 
p = 0.3397 e“/A^ [ 3 ], very close to pw = 0.332 e“/A^, so 
that |F/( 0 )| ~ 0 . 


4 Reflectivity results 

Figure |3] shows the reflectivity measurements (symbols) 
of 16 DMPC membranes on silicon substrates at four (out 
of nine measured) different osmotic pressures. The curves 
have been stacked vertically for clarity, with increasing 
pressure from 4 kPa (bottom) to 195 kPa (top). The con¬ 
tinuous lines are simulations based on the model described 
above with the corresponding electron density profiles shown 
in Figure I3A. The reflectivity spectra presented in Figure 
[3] were each scaled by the corresponding Fresnel reflectiv¬ 
ity. The error bars in the intensity at point n, Z\/„, are 
estimated considering Poissonian statistics (both for the 
raw reflectivity and for the offset scan). The error bar in 

the qz direction is taken as Aqn = Aq = 5 x 10“'^A , 

corresponding roughly to the symbol size, and is given by 
the estimated precision in sample alignment. 

7 ,From the electron density profiles one can see that 
the increase of periodicity d with decreasing pressure is 
mainly due to changes in the water layer thickness, while 
the bilayer structure is essentially invariant over the range 
in n studied, with a headgroup spacing (distance be¬ 
tween the two maxima in the electron density profile) 
duH = 34 ± 0.5 A, in good agreement with the value of 
34.4 A given by Petrache et al. [3]. The simulations match 
the measured reflectivity curves not only at the position 
of the Bragg peaks, but in the whole continuous 9 -range of 
the measurement. At lower osmotic pressure TT the higher 
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Fig. 3. Left : Symbols : measured reflectivity curves of 16 DMPC membranes on silicon substrates under different osmotic 
pressures, normalized by the corresponding Fresnel reflectivity. Solid lines : show full g-range simulations using the described 
model, with the density profiles given in Figure!?] The polymer concentration and the osmotic pressure are specified for each 
curve. Right : Residues of the reflectivity data, normalized by the estimated standard deviation (see text for details). The graph 
is horizontally split into two panels, with different y axes : Below 0.25 A the large discrepancies correspond to systematic 
errors, while for qz > 0.25 A the residues are more randomly distributed and their scaled amplitude is below 5. 


pgj electron density [eVA^] 



-30 -15 0 15 30 

z[A] 


Fig. 4. Electronic density profile of the bilayer, as a function 
of the applied osmotic pressure U. 


order peaks are suppressed due to increased thermal fluc¬ 
tuations, as quantified by the above model. As an illustra¬ 
tion, Fig. |S]4 shows the increase in the fluctuation ampli¬ 
tude as a function of the membrane index for a 16 bilayer 
sample. 

As discussed above, the relevant parameters for the 
reflectivity curves are the mean squared fluctuation am¬ 
plitudes (u^), which give access to the Caille parameter 
77. In order to compare our data to the bulk results [3] 
we then compute the interbilayer spacing fluctuation am¬ 


plitude cr^ = ((m„ —= rj—. Soft confinement 

TT^ 

theories predict an exponential dependence [3] of param¬ 
eter with the interbilayer distance, which can be taken 
as the thickness of the water layer, given hy dw = d — ds, 
where ds = 44 A is the thickness of the DMPC bilayer [3]. 

Fig. |SJ3 shows the reciprocal of the fluctuation ampli¬ 
tudes as a function of the water spacing cf~‘^{dyj), along 
with a fit to an exponential decay (solid line). The data 
points can be fitted to an exponential function oc 
exp(—du,/Afl), with a decay length Afl = 4.1 ± 0.9 A, com¬ 
parable to the results of the bulk study, which also re¬ 
ports exponential behaviour (shown as a dashed line) with 
Afl = 5.1 A [3]. 

The lamellar periodicity d was measured for all values 
of the osmotic pressure 7T, up to 870 kPa. We show n{d) 
in Figure | 6 ] (diamonds). For comparison, we also plotted 
the fit by Petrache et al. of the bulk data, Figure 7, up¬ 
per panel in their paper [3] (dashed line). They performed 
the measurements at 30 °C and obtain do = d{n = 0) = 
62.7 A, while our experiments, performed at 40 °C, yield 
do = 61.9 A. This deviation is in agreement with recent 
data of do{T) for DMPC [26]. A quick test can be per¬ 
formed by shifting their curve by 0.8 A to lower d values 
(solid line). The agreement is good, but a more detailed 
analysis is obviously needed for a meaningful comparison. 


5 Interaction potentials 

We saw above that our data can be brought to agreement 
with the bulk equation of state 77(d) by the 0.8 A shift, 
which can be attributed to the temperature difference. 
Consequently, the interaction potentials derived from the 
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Fig. 5. A) Amplitude of the bilayer position fluctuations (w^) in a solid-supported membrane stack consisting of 16 DMPC 
membranes at partial hydration, controlled by an osmotic pressure of 4 kPa (used for the model in Fig. [21 bottom). B) Reciprocal 
of the squared fluctuation amplitudes a~^{d) = tP 'as determined from the reflectivity fits, and exponential fit (solid 
line). Also shown is the exponential fit to the bulk data of Petrache et al. [3]. For clarity, their experimental data points are not 
shown. 



d[A] 


Fig. 6. Osmotic pressure 77 as a function of the lamellar d 
spacing. Diamonds : experimental data points. Dashed line : 
Fit to the bulk data of Petrache et al. [4]. For clarity, their 
experimental data points are not shown. Solid line: same curve 
shifted to lower d values to account for the effect of temperature 
(see text for details). 


bulk and the present data will also be identical or at least 
similar. However, the choice of the functional forms for 
the potentials, the geometric partitioning of the bilayer 
to calculate the van der Waals part and the choice of the 
bending rigidity k can all be debated. We therefore first 
give a brief discussion of the different interaction poten¬ 
tials used for neutral lipid membranes, and then present 
results based on modelling the equation of state n{d). 

The hydration potential is usually empirically described 
by an exponential function of the water layer thickness 

m ■■ 


(5) 


with a prefactor of the order of Hq = ksT /A^ and a decay 
length on the order of a few Angstroms A = 1 — 2 A. For 
the van der Waals potential, Petrache et al. [3], used the 
following expression and geometric convention: 


^vdw{diu') — 




vdW 


IGtt 


1 


1 


{dw + ds)^ (du, -|- 2dij)^ 

( 6 ) 

where d^, is the water layer and ds = d — dw the bilayer 
thickness. The expression should be regarded as an ap¬ 
proximation to the result of a more detailed treatment, as 
discussed by Fenzl m- Accordingly, the potential should 
be calculated from 


In 1 


Fvdw{dh, T) = 0.9 


keT 

Sirdl 



dxx 


- exp{-ax/dh)) \" . . 

l-Alexp{-ax/dh) ) ^ ^ 


(7) 


where dh is again the thickness of the hydrophilic lay¬ 
ers consisting of the water layer and the headgroups and 
An = (eH 2 o(w„) — ecH 2 (wn))/(eH 2 o(wra) + ecH 2 (wra)) is a 
function of the frequency-dependent dielectric constants 
of hydrocarbon and water. The prime symbol ' indicates 
that the static term (n = 0 ) has to be multiplied by 
1/2. The calculation is somewhat involved, however Fenzl 
has shown that a frequently used approximation of Eqs. 
[5] is valid for the dispersion term, but not for the static 
term which dominates under salt-free conditions. More¬ 
over, Podgornik and coworkers have shown that nonpair- 
wise additive contributions to the van der Waals interac¬ 
tion play a significant role in multilayers at large swelling 
[28j . However, for the present parameters, the above treat¬ 
ment should be sufficient. 

Apart from the molecular forces discussed above, steric 
forces resulting from membrane bending elasticity should 


fhydi.dw') — 77o exp( dw jF) 
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be included, as first introduced by Helfrich m- Accord¬ 
ingly, a repulsive undulation force arises 


fui = 0.42 


Kdl 


( 8 ) 


which cannot, however, be simply added to the molecular 
forces. Instead, steric forces must be treated by field theo¬ 
retical approaches which go beyond the mean field approx¬ 
imation [T], or by self-consistent models [30], but which 
to date have not been combined with realistic molecu¬ 
lar potentials in multilamellar stacks. Facing these com¬ 
plications, Petrache and coworkers [3] have pointed out 
that the measured rms-fluctuation of the next neighbor 
distance a = y^d/7r can be used to experimentally de¬ 
termine the fluctuation pressure Pfiuct, which they then 
added to the pressure calculated from the molecular po¬ 
tentials to fit their data n{d) = Pmoi + Pfiuct- Obviously, 
this approach avoids the problematic identification of the 
thermodynamic compression modulus to the bulk mod¬ 
ulus B as defined in the Caille model, but still assumes 
that the total pressure can be written as a sum, which 
may strictly only be true in mean field approximation. 
The advantage of the approach is that it makes clever use 
of the experimental information from either of the inter¬ 
connected functions B{d), r][d), or a{d) to compute the 
pressure. According to [3] 


(AksT)^ 1 da ^ 
(87r)2 K d {dw) 


(9) 


Fig. [5] shows the measured parameter the inverse of the 
fluctuation amplitudes a~^{d), along with a fit to an ex¬ 
ponential decay (solid line). The data points can be fit¬ 
ted to an exponential function 2.14 exp(—du,/4.17) A 
Subtracting the corresponding fluctuation pressure Pfiuct 
obtained by differentiation according to Eq. [3] for a given 
parameter n from 11 (d), the molecular interactions (hy¬ 
dration and van der Waals interactions) can then be mod¬ 
elled and compared to the data, as shown below in Fig. 
0 

Below we give results for two different approaches in 
the data analysis. The first approach is described in detail 
in |31) . It is based on the assumption that the periodicity d 
is dominated by the molecular forces, and that steric forces 
are comparatively small for relatively stiff phospholipid 
membranes. 

First approach: The calculation of the van der Waals 
interaction was based on equation © for the static con¬ 
tribution. The static part for n = 0 was numerically inte¬ 
grated between 0 and 100. For water eH2o(0) = 80 and 
for hydrocarbon (tetradecane) ecH 2 ( 0 ) = 2 was taken. 
For the dispersion term a hydrophobic bilayer thickness 
26 A and a Hamaker constant iJdis = 0.297 was used. The 
latter value has been chosen to approximate evalu¬ 
ated for dispersion relations eH 2 o(w) and ecH 2 ('^) which 
have been parameterized by oscillator models as in HU. 
Note that in this approach there is no free parameter for 
the vdW interaction. Figure [T] B shows the total interac¬ 
tion potential (solid line) with separately calculated static 


and dispersion terms as described above. A calculation 
based only on Eq. |6] (with adjusted Hamaker constant) is 
shown for comparison (dashed curve). The parameters of 
the hydration interaction which were then freely adjusted 
were Hq = 4.8 ksT/A^ and A = 1.88 A. The 10% reduc¬ 
tion from the fixed values in m in the van der Waals term 
for equation © may be attributed to the fact that Helfrich 
repulsive forces have not been included in the force bal¬ 
ance. At the same time, it is interesting to note the value 
obtained for k in this approach. To estimate k we note that 
the Caille parameter r] = 2 ^ /d been determined 

from the full g-range fits to the reflectivity curves at differ¬ 
ent pressures. As the compression modulus B{d) = dPL/dd 
has been independently determined by numerical deriva¬ 
tion of the measured 11 (d) curve, one can now estimate the 
bending modulus k from the experimental values of if^d). 
The best agreement was obtained for k = 23.2 ± 2.5 ksT. 
Within these uncertainties n ~ 2‘ikBT compares quite 
well with the value of k = IQfcsT at 30 °C determined 
from bulk suspensions of DMPC by Petrache et al. [3]. 
Note however that another study employing full g-range 
fits and osmotic pressure variation reports k = 11.5 fcsT 
m, again at full hydration and comparable T. Finally, 
thermal diffuse scattering analysis points to significantly 
smaller values k = 7 [33]. Note that the determina¬ 

tion of K from the osmotic pressure series is based on a 
problematic assumption, i.e. that the identification of the 
bulk modulus B as defined in the Caille model and the 
thermodynamic compression modulus is correct. More de¬ 
tails on the data analysis following this approach based on 
molecular interactions only are given in m- 


In the second approach we followed exactly the pro¬ 
cedure given by |3]. First, the fluctuation pressure is sub¬ 
tracted from P{d), and then the resulting bare pressure 
is modelled in terms of the molecular interactions. How¬ 
ever, the calculation of the fluctuation pressure according 
to equation|3]needs the bending rigidity k as an additional 
parameter. To illustrate the range of parameter variabil¬ 
ity, we present a comparison of two different choices of 
parameters: (a) all parameters and functions are kept as 
close as possible to those used in [3], in particular keep¬ 
ing K = 18.5 fixed. The corresponding values for 
Hq — 0.020 J/m^ and Hydw = 4.91 • 10“^^ J are practi¬ 
cally identical to the values in |3], showing that the same 
approach can explain both bulk and thin film data. The 
same treatment has then been carried out for a different 
choice of K = 8 fcsT. Again the simulations can be brought 
into agreement with the data, but only for a different set 
of parameters Hq — 0.028 J/m^ and Plvdw = 1.1 • 10“^° J. 
Thus values at the lower and upper range of the n values 
reported in the literature both lead to reasonable agree¬ 
ment, indicating that extra information from other ex¬ 
periments is needed to unambiguously determine the po¬ 
tentials. The potentials are shown for the two cases, and 
can be compared also to the potential in Fig. 1(b), de¬ 
rived from the data analysis under the assumption that 
the fluctuation repulsion is negligible. 
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Fig. 7 . Osmotic pressure 11 as a. function of lamellar peri¬ 
odicity d (same data as in Fig. EKA)). The simulations have 
been carried out following the approach of [3] for (a) fixed 
K = 18.5/cbT and (b) k = SksT. The parameters of the 
hydration interaction Ho (Eq. [Sj and the Hamaker constant 
Hvdw (Eq. O are varied to match the data. The simulations 
in (a) and (b) show the total pressure (solid line), the fluctua¬ 
tion pressure Pfiuct (dotted line), the pressure corresponding to 
the molecular potential (dashed line), and the potential trough 
corresponding to the molecular forces (dash-dotted line). The 
fluctuation pressure as determined from Eq.E]and the data in 
Fig. [5] was added to the pressure stemming from the molecular 
forces. 

6 Discussion and Conclusions 

In conclusion, we have presented an osmotic pressure ex¬ 
periment on thin solid-supported lipid multilayers (oligo- 
membranes). The x-ray reflectivity has been measured and 
modelled over the full ^z-range up to 0.7 A From this 
analysis fluctuation and structural parameters can be ob¬ 
tained, similar to the lineshape analysis of bulk suspen¬ 
sions [5115] . Solid-supported oligo-membranes offer some 
advantages both over thick multilamellar films and the 
bulk counterpart. Long range thermal fluctuations are not 
as strong as for bulk samples, and the scattering can be 
probed up to higher momentum transfer. Owing to the 


smaller number of bilayers, destructive interference in- 
between the Bragg peaks is not quite as strong as in 
thick stacks of several hundred bilayers. To achieve sat¬ 
isfactory hts, two important effects were taken into ac¬ 
count: (i) the static defects leading to a decreasing cov¬ 
erage of the bilayers from the substrate to the top of the 
him, and (ii) the thermal fluctuations of the bilayers sub¬ 
ject to the boundary condition of a flat substrate |24) . 
While (i) most likely reflects non-equilibrium aspects of 
sample deposition and/or equilibrium wetting properties 
(not further analyzed here), (ii) is exploited to deduce in¬ 
teraction parameters in the framework of linear smectic 
theory. The curve B(d) derived from the osmotic pressure 
series 77(d) is subsequently modelled based on different 
interaction potentials. However, this modelling cannot be 
carried out without assumptions or additional theoretic 
arguments. 

In the data analysis, we have presented two entirely 
different approaches to illustrate how the determination 
of interaction forces depends on the specific assumptions, 
theoretical arguments, or extra information taken from 
other experiments. The first approach builds upon the 
rather strong assumption that steric forces are negligible 
and that the derivative of the equation of state dU (d) /dd 
can be identified with the modulus B{d) which controls 
the thermal fluctuations. It then yields the parameters of 
the hydration force necessary to balance the van der Waals 
attraction at each given osmotic pressure. This approach 
also gives a value for the bending constant from simulta¬ 
neous inspection of B(d) and 77 (d). However, the resulting 
K ~ 23kBT, is probably an overestimation. The rather 
large value may point to the fact that B is underestimated 
by the contribution of only the bare potentials. Adding a 
fluctuation pressure would tend to increase B and thus 
decrease k = K d. Note that this determination of k is 
conceptually very different from a more direct assessment 
of K, e.g. from the measurement of diffuse scattering. 

The second approach includes the steric Helfrich un¬ 
dulation forces. This contribution is a subtle issue for the 
following reasons: (a) it has been shown by Lipowsky and 
coworkers that the Helfrich term cannot be simply added 
to the molecular forces. If one nevertheless uses a mean 
field approach, (b) the functional form to be used as well as 
the numerical prefactor are still under debate [34] . There¬ 
fore, we have followed an idea of Petrache et al. [3], who 
have carried out an osmotic pressure study on DMPC, 
which is the bulk analogue of the present work. Calculat¬ 
ing the partition function within the linear smectic elastic¬ 
ity model, they have derived derived an expression for the 
fluctuation pressure as a function of a measurable quan¬ 
tity, namely the derivative of the fluctuation amplitudes, 
see Eq. O In a mean field treatment, they add this pres¬ 
sure to the bare pressure calculated from the interaction 
potentials and fit the sum to the measured curve 77 (d). In 
this step, an assumption of k has to be made, e.g. from 
other experimental data. This approach has been carried 
out for two choices of k, see Fig. [T] 

We point out, however, that the questions related to 
the interaction potentials arise only on a secondary level. 
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where structural results {d and a) are interpreted and 
transformed to elasticity and interaction parameters. On 
the primary level that the structural results presented 
here, i.e. n{d), p{z,n) and (u^) (d) are well supported 
by the curves and fits shown here. The second level is 
necessarily model-dependent. We have presented two al¬ 
ternative approaches to illustrate the relation and interde¬ 
pendencies of different assumptions and results. It may be 
justified to conclude that the second approach as proposed 
and used by [3] is more appropriate, since steric repulsion 
is known to be important. However, the choice of the van 
der Waals expression may have to be improved to a more 
accurate form, and the choice of k is also an important 
issue. Unfortunately, the second approach also relies on 
the validity of a mean field approximation. This simplifi¬ 
cation could be eliminated in the future by generalization 
of a recent self-consistent calculation for bilayer fluctu¬ 
ations and interactions m to the case of several mem¬ 
branes or by use of the approach developed in [3S]. Fur¬ 
thermore, non-linear effects due to the asymmetry of the 
potential could also be included by more general models 
[5S1I55] and/or numeric simulations. To elucidate the valid¬ 
ity of the linearized model a posteriori, the rms-deviation 
\/{{un — Un+i)'^) between neighboring membranes can be 
compared to the width of the inter-bilayer potential well, 
see Figured] (B). For TV = 16 and 77 = 0.08 (full hydration) 
the bilayers in the center of the stack already exhibit con¬ 
siderable next-neighbor distance fluctuations in the range 
of 4 — 5 A, when compared to the water layer thickness. 
Thus the errors made in the simplifying assumptions are 
probably not negligible. We note, however, that at least 
under high osmotic pressure, where fluctuations are small, 
both the mean held approach and the harmonic approxi¬ 
mation for the potential should hold. 


Appendix 

In this appendix, we give a sketch of the derivation of 
our formula dm) for the reflectivity, insisting upon the 
separation between the specular and diffuse components. 
This is a classical result and a more detailed derivation 
can be found in references [37] (equation 2.28) and [33] 
(subsection 3.8.3) for the case of single interfaces and in 
[33] (section 3) for multiple interfaces. 

It is well known that bulk lamellar phases exhibit the 
Landau-Peierls instability, leading to a characteristic power- 
law variation of the scattered signal [5|. In such a system 
the fluctuation amplitude (m^) diverges. It is then more 
appropriate to use the correlation of the height difference, 
which remains finite for all finite values of r : 

5m»(r||) = ((wm(r||) - Wn(0))^^ (10) 

= (Wm) + (Wn) - 2(u™(r||)u„(0)) 

It is then easy to show [371[35] that the structure fac¬ 
tor of the lamellar stack (without taking into account the 


substrate contribution, so only the third term in Eq. (EU 
is described) reads : 

S'(q) = ^ ^-iq.d{m-n) f g-iq||r|| g- (r) ^ 
m,n 

where r = |r|| |, assuming that the fluctuations are isotropic 
in the membrane plane. For bulk systems, lim granii") = 

r—^oo 

00 , so that lim exp \—hq^gmn{'r)\ = 0 and the Fourier 

- - >00 ^ J _ 

transform with respect to ry in formula im yields a ” smooth’ 
function S'(q||) at hxed Pz. If, however, gmnir) does not di¬ 
verge for r —>■ 00 , the function exp [—^qlgmn{r)] now has 
a constant background, at a value of exp [— ^q^gmuioo)] = 
exp [—^ql ((Wm) + (^n))]i quantifying the ’’remanent or¬ 
der” in the system. Its Fourier transform is a Dirac delta 
function (5(q||) (in practice, its width is given by a com¬ 
bination of resolution effects, beam coherence and system 
size). This term is sometimes called the ’’true specular 
component”, because the smooth function discussed above 
(the ’’diffuse” component) also contributes to the specu¬ 
lar signal S'(q|| = 0, q^). However, as the diffuse scattering 
varies over a much larger qy range, it can be accounted 
for in the first approximation by an offset scan taken close 
enough to the specular sharp peak (see figure [3]). Finally, 
we can write : 

5'spec(q|| = 0,qz) = 

m,n 

( 12 ) 

which is the form employed in equation (13.11) . It is note¬ 
worthy that this ’’true specular” contribution is distinct 
in nature from the signal measured in SAXS experiments 
on powder samples, where only the diffuse signal persists. 
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